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Abstract. For a few nearby active nuclei, the disc orbiting the black 
hole is traced by water maser emission. By combining a simple model 
together with observed velocity profiles, we show that it is possible to 
put constraints on the black hole mass and on the distribution of matter 
(shape, density, size) in the outer disc. We then report possible parame- 
ters for the non-keplerian disc and for the black hole in NGC 1068, and 
mention an uncertainty of at least 25% on the central mass in NGC 4258. 



1. The inverse problem 

Finding the spatial distribution of matter p(f) that reproduces a given rotation 
curve n o bs. if) is a common inverse problem in Astrophysics. It can be addressed 
at the scale of AGN gaseous discs if one interprets the rotation curve of water 
masers detected in the core of some objects (Greenhill, 2002) in terms of a 
pure gravitational attraction. Efforts must therefore be made in developing 
reliable theoretical and/or numerical techniques to perform the inversion (e.g. 
Sibgatullin, Garcia & Manko, 2002). In the present situation, this must account 
for long-range interacting components (black hole, stellar cluster, bars, galactic 
bludge, etc.) as well as for objects populating the parsec-scale (disc, clouds, 
compact objects, torus, etc..) where masers are seen. Uncertainties are large, 
limiting our predictions and interpretations. We discuss here an example of 
inversion by considering a very simplified, two component system, made of a 
black hole and a disc, with an application to NGC 4258 and NGC 1068. 



2. Assumptions, equations and method 

At the parsec-scale, the standard theory predicts that all discs are (strongly) 
self-gravitating (Collin & Hure, 2000) and so, they might evolve towards a state 
with a unity Q-Toomre parameter. Accounting for disc mass effects is there- 
fore necessary to investigate these outer regions. But a reliable modeling of 
such discs is not permitted yet, given difficulties of various origins (nature of 
the flow, turbulent transport, numerical resolution, instabilities, boundary con- 
ditions), and some assumptions are necessary. The simplest way to tackle the 
question is a "minimal approach" which assumes that the velocity field has a 
pure gravitational origin. For a black hole with mass Mbh and a disc only, the 
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most basic equation set to solve is 
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where g is the gravitational acceleration due to the disc, R is the cylin- 
drical radius, r is the spherical radius and p is the unknown density, a tri- 
dimensional function. In other words, if we are able to compute the self- 
gravitating field for discs with various size and mass distributions, we can select 
among possible solutions those matching observations. Note that, in this kind 
of approach, no physical model for the disc imposes the p— function. 

Here, the Poisson equation is solved numerically assuming azimuthal invari- 
ance for a given p{r) from the integral form 



by the method described in Hure (2002a). We stress that the full regularization 
of the volume integral requires three successive integrations, otherwise diverging 
fields (and subsequently infinite velocities) can result at the disc edges (e.g. 
Mestel 1963). This can be critical if the maser spots rotate near the disc outer 
edge. The present inversion method works with three inputs: the black hole 
mass Mbh, the mass density in the disc p(R, z), and its outer radius R ou t- As 
several solutions can lead, at least locally, to similar rotation curves, we impose 
that the total surface density £ = j pdz is of the form S oc R s , with a smooth 
connection to the ambient medium at the disc edges (see Hure 2002b for details). 
The mass density inside the disc obeys a parabolic law with the altitude, and the 
disc aspect ratio H/R is held constant. These additional restrictions however 
do not remove the degeneracies. Once a model fits observations, inputs become 
outputs. 

3. NGC 4258: the possible ambiguity 

NGC 4258 is the most famous example of Seyfert-2 active galaxy showing a thin 
outer disc with masers in quasi-perfect keplerian motion (Myoshi et al. 1995). A 
priori, this sustains the simple idea that most of the mass is effectively gathered 
at the center, implying a black hole mass ~ 3.9 x 10 M Q . If this is correct, the 
inversion method can furnish an upper limit for the disc mass: as figure 1 shows, 
the mass of the masing disc at ~ 0.28 pc (or ~ 8 mas) must be less than ~ 10 5 
M© (that is a mass ratio q < 2.5 x 10~ 3 ), otherwise a deviation to the keplerian 
rotation results. But the possibility that a massive disc is present around a 
lighter black hole must not be ruled out. Actually, the rotation curve concerns 
a short spatial range and various i?-laws can mimic a keplerian profile. Figure 
lb shows a case where the rotation curve is still well reproduced with a black 
hole 25% less massive and a disc with a mass ~8x 10 6 M (or q ~ 0.26) and 
with an outer edge at 0.11 pc (in this case, the masing clouds rotate "outside" 
the disc). In this example, we have deliberately selected a disc with Q ~ est (for 
the reason mentioned before), but the upper limit derived in this way is almost 
not affected by this specific choice. 
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Figure 1. Velocity of masers versus the angular separation in NGC 
4258. Left: theoretical velocity expected for a black hole with nominal 
mass M B h ~ 3.9 x 10 7 M and discs with mass of 10 7 , 10 6 and 10 5 M© 
and outer edge at ~ 8 mas. Right: Same but with a 25% lighter black 
hole and with a disc with mass ~ 8 x 10 6 M (see text and Tab. 1). 



4. NGC 1068 

In contrast, the rotation curve of the outer disc in NGC 1068 appears in no- 
ticeable departure with respect to the keplerian case. Greenhill et al. (1996) 
found ^obs. °c -R -1 ' 31 in the range 0.65 — 1.1 pc, and mentioned that the disc 
self-gravity could be the possible source of this non-keplerian behavior. Other 
mechanisms can be invoked, like radiation pressure effects (Pier & Krolik 1992) 
or gravitational attraction by a stellar cluster (Kumar 1999), but the disc self- 
gravity appears as a natural explanation. Applying the inversion method gives 
several families of solutions. The sensitivity to s is important. Two possible so- 
lutions characterized by a mean angular rotation index of ~ —1.31 are shown in 
Fig. 2. The parameters deduced for the disc and central black hole are reported 
in Table 1. We find that the disc mass is > 10 7 M , comparable to that of the 
black hole. We stress that surface densities that do not vary approximately as 
Rr l do not match observations: with s > — 1, the global rotation curve is flatter 
(and steeper with s < —1). The important point is that the value inferred for 
the Q— parameter is of the order of unity, meaning that the outer disc could be 
in a marginally stable state. Moreover, the temperature and the density in the 
outer disc seem appropriate for excitation of water molecules (Hure 2002b). 



parameters 


NGC 4258 
light disc massive disc 


NGC 1068 
thick disc thin disc 


M bh , in M 

i?out ; in pc (in mas) 

disc aspect ratio H/R 

index s 

q(Rout) 


4 x 10 Y 3 x 10 Y 
~ 0.28 (8) > 0.11 (3.1) 
0.003 ' 0.1 
-1 -1 
< 0.0025 0.26 


1.2 x 10 Y 1.2 x 10 Y 
> 1.3 (18) > 1.5 (21) 

0.3 «0.3 
~ -1.05 1 

0.91 0.65 



Table 1. Possible parameters for the disc and for the central black 
hole in NGC 4258 and NGC 1068 found from the inversion method. 
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observations 
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Figure 2. Same as for figure 1, but for NGC 1068. Left: a thick disc 
solution. Right: a thin disc solution (see Tab. 1). 



5. Conclusions 

We have discussed possible constraints (and uncertainties) on the mass of the 
central black hole and on structure of the parsec-scale region in NGC 4258 and 
NGC 1068, by a method involving a simple axi-symmetrical disc model and a 
black hole. The subsequent challenges will then be i) to understand, in terms of 
physical processes, why the disc has the inferred properties, ii) to check the the 
minimal approach, and hi) to discuss the implications of these solutions for the 
inner disc (see Lodato & Bertin in these proceedings). 
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